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Axisymmetric force-free magnetosphere of a pulsar. I. The 
structure close to the magnetic axis 
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ABSTRACT 

The stationary axisymmetric force-free magnetosphere of a pulsar is studied analyti- 
cally. The pulsar equation is solved in the region close to the magnetic axis. Proceeding 
from linearization of the current function in the axial region, we find the axial mag- 
netic flux function valid at any altitude above the neutron star. This function is used 
as a starting approximation to develop series for the non-linear pulsar equation in 
the polar region. Taking into account the quasi-monopolar character of the pulsar 
magnetic flux at infinity, we obtain unique asymptotic series for the flux and current 
functions. At infinity, both functions are close but not equivalent to those known for 
the case of a force-free monopole. The flux function at the top of the polar gap is 
found to differ from the dipolar one at the neutron star surface. With our results, the 
transverse current sheet closing the pulsar circuit at the neutron star surface is con- 
sistently incorporated into the global magnetospheric structure, the backward particle 
flow at small polar angles can be excluded and the stationary cascade scenario looks 
admissible. The present paper is the first step toward complete analytic description 
of the pulsar force-free magnetosphere allowing for the plasma-producing gaps and 
pulsar current circuit closure. 
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1 INTRODUCTION 



^ ■ Pulsar magnetosphere 



contains 
which 



the plasma 
sufficiently abun- 



l|Goldreich fe Julian! Il969l) . wnicn is 

dant to affect the magnetic field structure. This necessitates 
a self-consistent consideration of fields and currents in 
pulsars. The basic model of the pulsar magnetosphere is 
that of a rotating axisymmetric force-free dipole, where 
the magnetic and rotational axes are aligned, the elec- 
tromagnetic forces are balanced and the particle inertia 
is ignored. Then the poloidal current and magne tic flux 
are re l ated by the well-known pulsar equation ((Michel 
Il973al : IScharlemann fe Wagoner! 1 19731 : lOkamotol 1 1974 ). 
Being formulated almost 40 years ago, it still lacks a proper 
analytic solution for the dipolar case. 

An ex act solution is known only for the force-free 
monopole (Michel 1973 aj), including the vers i ons o f a split 
and offset monopole <|Miche]| fl99lT iPetroval l2012t ). which 
are believed to replicate the behaviour of the force-free 
dipole at large distances from the origin. In the case of 
a monopole, the problem is strongly simplified, since the 
poloidal component of the force-free magnetic field appears 
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to be the same as in the original vacuum configurations. 
In the case of a dipole, the flux and current functions 
entering the pulsar equation are both unknown, and the 
problem is to guess such a current distribution that makes 
the self-consistent flux function obey the physically mean- 
ingful boundary conditions. The early attempts of solving 
the pulsar equation for the current functions of a special 
form yielded the flux functions which were valid only in- 
side the light c ylinder and could not be smoothly contin- 
ued t o infinity (Michel 1973b]; iBeskin. Gurevich fc Istominl 
1983). Later on a more sophisticated current function pro- 
viding a smooth passage of the magnetic field lines across 
the light cylinder was found by means of numerical meth- 
ods (|Contopoulos. Kazanas fc Fendtll 19991 ). The necessity to 
revise the set of commonly u sed boundary conditions was 
pointed out in IPetroval (|2012f ) . 

The pioneering work of IContopoulos et all (|1999T ) 
inspired extensive numerical studies of the pulsar 
equation. The numerical solution for the stationary 
axisymmetric force-free magnetosphe re was confirmed 
by the time-de p endent simulatio ns l|Komissarovl 120061 : 
iMcKinnevI 120061 ; ISpitkovskvl | 2006) and generalised in a 
number of aspects. In particular, the non-axisymmetric 
3D magnetosphere was considered (|Spitkovskvl |200d : 
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Kalapotharakos fc Contopoulosl 120091; Bai fe Spitkovskvl 
2010l ; Kalapotharakos. Contopoulos fc Kazanasl 20121 ). 
the differential rotation was taken into account 



(1 



dContopoulos! 120051; iTimokhin 20061. 20071), the non- 


ideal MHD case was addressed 1 


Kalapotharakos et al.ll2012l; 


iLi. Spitkovskv & Tchckhovskov 


120121). However, the set 



of boundary conditions used in all these simulations did 
not allow for the presence of the plasma producing gaps. 
Therefore the numerical treatment of the pulsar equation 
was not self-consistent. 

In our previous work l|Petrovall2012h , we have suggested 
an empirical model of the pulsar axisymmetric force-free 
magnetosphere, which for the first time includes the polar, 
outer and slot gaps into the global magnetospheric struc- 
ture. The basic prescriptions of the model are as follows. 
Although the polar gap height can be regarded as infinites- 
imal, the magnetic field structure at its top differs from the 
original dipolar one at the neutron star surface. As the outer 
gap is entirely located within the light cylinder, beyond the 
light cylinder the null line should not intersect the magnetic 
field lines. The polar and outer gaps control different bun- 
dles of the open field lines and produce the direct and return 
currents, respectively. Therefore beyond the light cylinder 
the null line coincides with the critical field line, both going 
parallel to the magnetic equator at a certain altitude above 
it. The compensation current necessary to close the pulsar 
circuit originates in the slot gap located along the critical 
field line inside the light cylinder. 

Our qualitative model needs to be quantified. With the 
present paper, we begin a systematic analytic description 
of the pulsar magnetosphere within the framework of this 
model. The present study is focused on the region close to 
the magnetic axis. The absence of the poloidal current along 
the magnetic axis provides an essential simplification: in the 
nearest vicinity of the axis, the generally non-linear cur- 
rent function and the pulsar equation on the whole can be 
linearised. In Sect. 2, we solve the simplified pulsar equa- 
tion in the axial region and find the magnetic flux func- 
tion at any altitude z above the neutron star. In Sect. 3, 
the complete non-linear pulsar equation is considered. We 
search for the solution in terms of series in p 2 /z 2 (where 
p is the separation from the axis) and use the axial flux 
function as a starting approximation. Taking into account 
the quasi-monopolar character of the flux function at infin- 
ity, we find a unique current function and the corresponding 
unique asymptotic form of the flux function at large alti- 
tudes above the neutron star. With the current function in 
hand, we also obtain the flux function at the top of the polar 
gap and demonstrate its distinction from that of a dipole. 
Our results and their implications are discussed in Sect. 4 
and briefly summarised in Sect. 5. 



2 SIMPLIFIED PULSAR EQUATION IN THE 
AXIAL REGION 

Our consideration concerns the model of a stationary ax- 
isymmetric force-free dipole. It is convenient to choose the 
cylindrical coordinate system (p, <j>, z) with th e axis along the 
pulsar axis. Then the pulsar equation reads ([Michel! Il973al : 
IScharlemann fc Wagonerll 19731 ; IOkamotolll974l )~ 



d 2 f 19/ SPj_ 

dp 2 p dp dz 2 



p dp 



(X) 



where the dimensionless functions f(p, z) and A(f) are pro- 
portional respectively to the magnetic flux and poloidal 
electric current through the circle of a radius p centered 
at the magnetic axis at an altitude z. As we are inter- 
ested in the region close to the magnetic axis, the rele- 
vant boundary conditions are as follows. Along this axis, 
the flux / and current A are both zero. At infinity, the 
field lines become radial, / = f(p/z), and resemble those 
of a magnetic monopole. According to th e current models 
(e.g IShibata, Mivazaki fc Takaharal |2002| ; lYuki fc Shibatal 
l2012h . the pulsar current circuit closes through the neutron 
star. Then out of the circuit, well below the stellar surface, 
the poloidal current is zero, the magnetic field is dipolar, 
/ = p 2 /(p 2 + z 2 ) 3//2 , and satisfies the multipolar equation 



d 2 f 



d 2 f 



dp 2 p dp dz 2 



(2) 



In the region of interest, p <C 1, we have / fa p 2 /z 3 . Just 
above the top of the transverse current sheet, the current is 
aligned with the poloidal magnetic field and remains con- 
stant along a fixed field line. Note that the jump of the 
poloidal current from zero to the value dictated by the global 
force-free configuration, which occurs at the neutron star 
surface bearing the transverse current, is compatible with 
the continuity equation, V ■ J = 0. 

The presence of the infinitesimally narrow polar gap just 
above the transverse current sheet nothing affects the for- 
mulation of the force-free problem. Indeed, both unknowns, 
A and /, are not changed in the course of the particle accel- 
eration, and the role of the gap is only to supply a sufficient 
amount of plasma necessary to sustain the force-free config- 
uration. Therefore in our consideration the top of the polar 
gap and, correspondingly, the bottom of the force-free region 
can be assumed to coincide with the top of the transverse 
current sheet at the stellar surface. In the present paper, 
we also neglect the potential drop in the polar gap and the 
resultant differential rotation of the magnetosphere. 

The current function on the right-hand side of Eq. |T} 
is unknown. It is believed to be non-linear, in order to pro- 
vide the current circuit closure, and close enough to that of 
a force-free monopole, AdA/df = 2/(1 - ,f//o)(2 - f/f ), 
where /o is a constant (see, e.g.. IContopou los 2005). Note 
that the solution of the pulsar equation with an a ssumption 
of a linear current function (Bcs kin et al.lfl983T ) appeared 
non-physical, since it could not continuously cross the light 
cylinder. Indeed, the solution of a non-linear problem in 
general differs essentially from that of a linear one. At the 
same time, in the nearest vicinity of the magnetic axis, the 
non-linear current function can be expanded into the Tay- 
lor series in /. Keeping in mind its similarity to the current 
function of a force-free monopole, one can expect that the 
linear term of the expansion is non-zero. Thus, close to the 
magnetic axis, the current function can be linearised. 

With the above assumptions, the pulsar equation is sub- 
stantially simplified and takes the form 



8^1 
dp 2 



_^ _ 19/ , 9V _ _ , , 

»« 2 P dp + dz 2 aJX[z 



zo) 



(3) 
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where a is an arbitrary constant, x( z — 2 o) is the step func- 
tion, zo stands for the polar gap altitude, and the terms ~ p 2 
are ignored. Strictly speaking, for these terms to be ignored 
the step function x( z ~ z o) should actually be somewhat 
smeared at a scale h 3> p. Thus, assuming a narrow trans- 
verse current sheet, h <C zo, we are restricted to consider 
only the nearest vicinity of the magnetic axis, p <^ih, which 
we call the axial region. 

Eq. allows separation of variables, in which case the 
solution can be presented as 



/ = 



C{\)v{p)u{z)d\, 



where the functions v(p) and u(z) obey the equations 
dv 



ldv 2 

XT _ ~T~ + A v 
dp 2 p dp 



0, 



du 
dz 2 



A u — —aux(z — zo), 



(4) 



(5) 



(6) 



A is independent of p and z, and C(A) is an arbitrary func- 
tion. Since at the magnetic axis / = 0, we have v(p) = 
pJi(Ap). In order to provide the dipolar character of the 
magnetic field at z < zo, f = p 2 /{p 2 + z 2 ) 3 ^ 2 , it is nec- 
essary to choose u(z) = exp(— Xz) and C(A) = A. In the 
Appendix we solve Eq. ^ for the smoothed step func- 
tion and show that well above the top of the polar gap 
u(z) — exp(— z\[)? 



/ 



P - J vVT^e 



a). Then the flux function reads 
"V/i, (7) 



where p 2 = A 2 — a and it is taken that Ji(Ap) w Ap/2. 

Eq. l [ 71) can be ca lculated using the integral 



l|Gradshtevn fc Rvzhildll98ch 



(a; + u ) 



a^-i--i« da . : 



/ 2u 

2 \7 



v-1/2 



(8) 



where T(y) is the Euler gamma-function, H v _ 1 / 2 (up) the 
Struve function, and Y v _ 1 / 2 (up) the Bessel function of the 
second kind. Differentiating both sides of Eq. JSJ with re- 
spect to p and making use of the recurrent relations for the 
Struve and Bessel functions yield 



with 



4> 



v 3/2 



+ ^[h 2 (vW)-r 2 (^)] 



(9) 



Thus, we have found the axial flux function of the rotat- 
ing axisymmetric force-free dipole. In order to confirm the 
validity of our technique of simplifying and solving the pul- 
sar equation, one can analogously apply it to the force-free 
monopole, in which case the exact solution is kno wn. The 
monopolar condition at the origin, / = 1 — z/y^z 2 + p 2 , 
implies that O(A) = 1. Then, instead of Eq. lff)l. we have 
2 



3 2 



z pdp, 



which is exactly reduced to / = p 2 /2z 2 , just as could be 
expected. 

Taking into account the approximations at z — s> 0, 



r(3/2)r(i/ + 3/2) ' 



: rMK/2)- 



one can see that at z < 1 the magnetic field structure is 
roughly dipolar, / ~ p 2 /z 3 (recall that we consider the re- 
gion very close to the magnetic axis, p/z <C 1). With the 
asymptotic relation 



m — 1 

H v (0 - Y„(0 = i V 

7T * — ' 



r(fc + i/2) 



T(v + l/2 - fc)(£/2) 2 *-"+i 



+ 0(\z\ 



(10) 



we find that at large altitudes, z ^> 1, the flux function bears 
the monopolar character, / ~ *Jap 2 /2z 2 . This is just what 
is expected for the force-free dipolar field and attributed to 
the action of the current, which flows through the polar gap 
and further along the open magnetic lines. Our consideration 
shows that, at least in the axial region, the flux function 
of a force-fre e dipole asymptotically coincides with f m on = 

y/Zi(l-z/y/F+*). 

The function <&(z) along with its asymptotes is plotted 
in Fig.[T] Figure [2] shows the magnetic field lines in the axial 
region for the cases of a force-free and vacuum dipole. One 
can see that the force-free field lines tend to go closer to 
the axis. This seems to be in line with the general fact that 
the force-free dipole has an additional bundl e of open field 
lines as compared to the vacuum dipole (e.g ., lMichellll973bl ; 
IContopoulos et al.|[l999l ; ,Contopoulosl,2005 



3 COMPLETE PULSAR EQUATION AT 
SMALL POLAR ANGLES 

In the previous section, we have obtained the axial flux func- 
tion for the rotating axisymmetric force-free dipole. With 
this result in hand, we shall treat the complete pulsar equa- 
tion |l} at small polar angles, 6 — atan(p/z) <C 1. The 
solution can be searched for in the form 



E 



, / \ 2k / 2k 
b k (z)p /Z , 



(11) 



where 61(2) = z 2 $(z). Substituting Eq. (iTTj) into Eq. |T} 
leads to the recurrent relation 



4fc(fc + l) 



9fc+i 

,2fc+2 



4k 



z 2k ^ z 2 \ z 2k J ^ z 2 y z 2k-2 j 



-AA'oi 



(12) 



where 



AA' 2k = —s 



k-l 



fe — 2 fc — i— 1 

EE 



bibjbk- 



(13) 



and £j are constants. The latter relation results from ex- 
panding A A' into the Taylor series in / and using Eq. (|lip . 

If the current function were given and, correspondingly, 
the set of £i were known, one would successively obtain the 
functions bk from Eq. (|12p . However, in our case the flux and 
current functions are both unknown and should be found 
self-consistently. For this purpose we additionally impose the 
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boundary condition at infinity. Namely, at infinity the mag- 
netic field lines should become radial, i.e. the flux function 
should only be the function of the polar angle, / = f(p/z), 
and, correspondingly, for any k 



lim bfe(z) = Constfc. 



(14) 



Applying this condition to Eq. (|12[) at each step of the it- 
eration procedure and using the asymptotic representation 
()10[) . one can obtain a unique value of £fc and the corre- 
sponding function bk+i(z), which is also unique. In practice, 
it is sufficient to proceed directly from the asymptotic form 
of h(z), 

i / \ _ 3 15 315 , _ _,, 

bl(2) " 1+ 4^-16^ + 64^---" (15) 
retaining fe terms in order to find the asymptotic value of 
bk+i- As a result, the current function of an axisymmetric 
force-free dipole and the asymptotic form of the correspond- 
ing flux function at z — > oo and p/z <C 1 read 

AA' = 4f-3f- 7 -f + ^-f i + ..., 

p 2 3p 4 lp 6 13p 8 9567 p 10 
1 z 2 4 z 4 2 z 6 32 z 8 1280 z 10 



+ . 



(16) 



Taking the axial flux function for the force-free 
monopole, / moB = p 2 /z 2 (in which case a — 4), we have 
exactly b\ — 1 and obtain analogously 



AA' mon = 4/ - 3/ 2 + -/ 3 , 
_ P 2 3p* 5f 

/mon — -o — 7— J + 



35/ _63_p^ 

64 z 8 + 128 z 10 + 



(17) 



z- 4 z 4 8 z 6 
One can see that in these two cases the flux and current 
functions differ only slightly, 

„6 



/ /mon 



*6 ' 



Ayr - AA 



-4/° 



(18) 



but far enough from the magnetic axis this distinction may 
well become substantial. 

We did not manage to guess the analytic form of the 
functions corresponding to the asymptotic series (fT6)) . While 
trying different transformations, we have only found that / 
can be presented in a more suitable form, 

/ = sin 2 + 7 sin 4 6>+|^ sin 8 + ..., (19) 



where 6 = atan(p/z). Since the term ~ sin 6 6 is zero, it is 
reasonable to keep only the first two terms of the series. 
Then at sin# = 1 we have / = 1.25, which is, unexpectedly, 
in a good corresponde nce with the value / = 1.23 obtained 
bv lContopouiosl (|2005h numerically. 

At large radial distances, r = yj z 2 + p' 2 — > 00, the 
asymptotic form of th e pulsar equation ^ reads (e.g., 
IContopoulos et al.lfl999h 

t 2 (l + t 2 )^ + (2t 3 +t)^=AA', 



(20) 



At 2 v 'At 
where t — p/z. Integration of this equation yields 

fVl + i 2 ^ = A. (21) 

One can check that the functions AA' and / given by 
Eq. (fTfJ)) do satisfy the relation (f2Tj) . 



Although the expression for the flux function obtained 
above is valid only at z — > 00, the expression for the current 
function is valid all over the magnetic field lines close enough 
to the axis. With this function in hand, from Eq. (|12l) one 
can obtain the flux function at z — > 0, i.e. directly at the 
top of the polar gap, 



/ 



_3ff_ 15 p b 203 p 8 

z 3 2 z 5 + 8 z 7 96 z 9 



+ . 



(22) 



(23) 



which differs from the dipole expansion as 

f- f Li 

J /dip gg z9 . 

By means of Eq. (|22|) . (the polar region of) the transverse 
current sheet at the neutron star surface is included into the 
global structure of the pulsar force-free magnetosphere. 



4 DISCUSSION 

Our paper presents one of the scarce attempts of analytic 
treatment of the pulsar equation. An idea to start such a 
consideration is motivated by the actual necessity to comple- 
ment the recent multitudinous numerical results on the pul- 
sar force-free magnetosphere with the physi cally grounde d 
theoretical advances. In our previous paper (|Petrovall2Qjl ). 
we have suggested an empirical model of the pulsar axisym- 
metric force-free magnetosphere, which includes the polar, 
outer and slot gaps. The present paper starts the analytic 
description of this model and concentrates on the region 
close to the magnetic axis. 

Since along the axis the current function is zero, in the 
closest vicinity of the axis it can be linearized. Consequently, 
in the axial region, the pulsar equation is essentially simpli- 
fied and can be solved by separating variables. As a result, 
we have obtained the axial flux function ((9]) valid at any 
altitude above the neutron star. It links the magnetic field 
below the stellar surface to that at the top of the polar gap 
and that at infinity. It is shown that at infinity the axial 
field of a force-free dipole asymptotical ly coinc ides with the 
field of a monopole, / mon = 2(1 — z/-\f z 2 + p 2 ) w p 2 /z 2 . 

The axial flux function ((9J is used as a starting approx- 
imation in the treatment of the complete non-linear pulsar 
equation in terms of series in p /z . Taking into account 
the quasi-monopolar character of the magnetic flux at large 
radial distances yields unique asymptotic series for the flux 
and current functions in the polar region. At z — > 00 both 
functions appear to differ slightly from those for the force- 



free monopole, / — f a 



-4f 



At the top of the polar gap, the flux function slightly differs 
from that of a pure dipole, / - /dip ~ 7p 8 /96z 9 . Although 
the asymptotic series obtained do not seem to be suitable 
for direct quantitative applications, our results have several 
implications. 

With the axial flux function in the form © , we have the 
second derivative in p at the axis, (d 2 f /dp 2 ) p =o = 2<I>(z), 
which can be regarded as an additional condition at a seg- 
ment of the boundary in the elliptic-type problem associ- 
ated with the pulsar equation. Within the framework of 
the quasi-reversal method of solving differential equations 
in partial derivatives, this may be used to compensate for 
the unknown part of the boundary or the unknown bound- 
ary condition. This is especially actual, since the full set of 
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boundary conditions in the pulsar force-free magnetosphere 
with the plasma-producing gaps is still unknown. 

The flux function at the top of the polar gap ()22p 
demonstrates that the magnetic field differs from the original 
dipolar one. Our result implies that it is the transverse cur- 
rent closing the pulsar circuit that is responsible for the mag- 
netic field distortion. Then an exact form of the transverse 
current is sufficient to successively find the magnetic field at 
the top of the polar gap, the force-free current distribution 
across the polar-gap-controlled field lines and the force-free 
field up to infinity. Thus, the long-standing problem of the 
polar gap adjustment to the global magnetospheric struc- 
ture is reduced to the problem of consistent incorporating 
the pulsar current circuit closure. 

Generally speaking, the global structure of the force- free 
fields and currents is a basis for examining the microphysics 
of particle flows all over the magnetosphere. Of course, such 
a detailed study is beyond the framework of the present pa- 
per, but with the flux and current functions obtained we can 
find the net charge and poloidal current densities which sup- 
ply the global magnetospheric configuration in the regions 
of interest. Retrieving the dimensional coefficients, we have 



B R 



0-f%S D 2 



2RI 



R 2 $(Z/Rl), AA' 



(24) 



where Bo is the magnetic field strength at the neutron star 
surface, Rns the stellar radius, Rl — c/Q the light cylinder 
radius, the angular velocity of pulsar rotation, R = pRl, 
Z = zi?L- Then for the charge and current densities, 



Pe 



n ( d 2 f i_df_ aV 

4ttc \dR 2 + RdR + dZ 2 



4-kRV VdR) \dZ) 



(25) 



we have 



Pe 



" 2ttc Rl 

_ QBo i?Ns 
Jp ~~2?rc Rl 



<t>(Z/R L 



$(Z/R L )c, 



(26) 



where — Q.Ba/2-wc corresponds to the Goldreich- Julian 
charge density at the top of the polar gap, the counteraligned 
rotator is chosen, it is taken into account that A' — —2/Rl 
and only the eldest terms in R are kept. As can be seen from 
Eq. (J26J, J p = p e c. 

The same result, J p /p e c — 1, is also obtained with the 
asymptotic flux function at infinity (|16[) . Note that, despite 
the distinction from the field of a monopole, the equality is 
fulfilled exactly, just as in the monopolar case. For the flux 
function at the top of the polar gap (1221) we have J p /p e c~ 
1 — p 2 /2z' i . As could be expect ed, this coincides with the 
relation (A2) in lTimokhinl ||2006l ) for the purely dipolar field, 
since the difference of the flux functions is very small (cf. 
Eq. (H3). 

Note that in the polar region considered in the present 
paper the poloidal current density differs only slightly from 
the Goldreich- Julian current density, < f — J P /p e c <C 1. 
This particular result is in line with the stationary pair cas- 
cade inside the polar gap and is likely to testify against the 
backward particle flow along the axis. Although the back- 
ward particle flow is thought of as a necessary attribute 



of the pulsar magnetosphere |Lvubarskvll2009l ), it rather 
streams along the field lines controlled by the outer gap. 
The recent numerical simulations of the coexisting gaps in 
the pulsar magnetosphere s how the backward par ticle flow 
along the axis (see Fig.l in lYuki fc Shibatall2012h . but the 
authors themselves admit that this, together with the non- 
zero accelerating electric field at the axis, is rather an arti- 
fact of their numerical method. 

Our consideration have demonstrated that the asymp- 
totic representations of the flux and current functions close 
to the magnetic axis are unique. In the polar region, these 
functions appear close eno ugh to those obtained by means 
of numerical simulations in lContopoulos et alJ l| 19991 ) and a 
number of subsequent works. This could be expected, since 
we have assumed the same boundary condition at the axis. 
Note that at larger polar angles the presence of the outer 
and slot gaps should markedly change the force-free field as 
compared to that known from the previous numerical simu- 
lations of the pulsar force-free magnetosphere. This will be 
a subject of our future study. 



5 CONCLUSIONS 

We have considered the stationary axisymmetric force-free 
magnetosphere of a pulsar and performed the analytic treat- 
ment of the pulsar equation close to the magnetic axis. 
Linearization of the current function in the axial region 
leads to the self-consistent flux function with the approx- 
imately dipolar and monopolar behaviour at low and high 
altitudes, respectively. The axial flux function has been used 
to construct the solution of the non-linear pulsar equation 
at small polar angles in terms of series. Keeping in mind 
the quasi-monopolar character of the flux function at in- 
finity, we have obtained unique asymptotic series for the 
flux and current functions, which differ from the monopolar 

ones as / - / mon p 6 /8z a , AA' - AA' mw 4/ 3 . With 

the current function obtained, the flux function at the top 
of the polar gap appears to differ from the dipolar one as 
/ — /dip ~ 7p s /96z 9 . This is attributed to the action of the 
transverse current closing the current circuit. Thus, our re- 
sult includes the polar part of the transverse current sheet 
at the neutron star surface into the global magnetospheric 
structure. 

Possible implications of our results can be summarized 
as follows. 

(i) The axial flux function found implies an additional 
boundary condition at the axis, (d 2 f /dp 2 ) p =a = 2<b(z). This 
may, at least partially, compensate for the unknown bound- 
ary segments and boundary conditions in the pulsar force- 
free problem. 

(ii) The self-consistent magnetic field at the top of the 
polar gap may constrain the physics of current closure in 
pulsars and give new insights into the polar gap theory. 

(iii) Our results testify against the backward particle flow 
at small polar angles. 

(iv) The relation of the resultant net charge and poloidal 
current densities does not exclude the stationary cascade 
scenario in the polar region of the polar gap. 

The results of the present paper are the first step in 
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the analytic descript ion of our mo del of the pulsar force- 
free magnetosphere which consistently in- 
cludes the plasma-producing gaps. As our model assumes 
the same axial boundary conditions as the classical one, in 
the polar region considered here our results are roughly in 
line with the numerical results obtained within the frame- 
work of the classical model. However, in other parts of the 
magnetosphere our model is essentially distinct, and the an- 
alytic description of fields and currents in these regions is of 
particular interest. This will be a subject of the subsequent 
papers. 
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APPENDIX A: CONTINUOUS TRANSITION 
ACROSS THE POLAR GAP 

We are going to find the solution of Eq. (JS]) at z 3> zo, taking 
into account that at z < zo u — exp(— Xz). It is convenient 
to introduce a smoothed step function in the form 



X(z - zo) = 



(z - z )/h 
l + (z-Zo)/h' 



(Al) 



Note that the results of this Appendix are not affected by 
the concrete way of smoothing the step function. The form 
(| All) is chosen for the sake of simplicity and it is not expected 
to have anything in common with the true distribution of 
current inside the transverse current sheet. With Eq. (|A1[1 
we reduce Eq. © to 



d 2 u 
dy 2 



(l + 2/)^i-^ 2 [A 2 + (A 2 -a) y]u = 0, 



(A2) 



where y = (z — za)/h. The solution of Eq. (IA2[) reads 



2h V / X 2 -a(y + lje-'v^" 
2h 



x ^ CiM 



1 + 



1 + 



2h 



2,2hy/X 2 -a(y + l) 



Vx 2 



,2,2h V / X 2 -a(y + l) 



(A3) 



where M(a, b, x) and U(a, b, x) are the degenerate hypergeo- 
metric functions. At h — > one can make use of the relations 



M(v + 1/2, 2v + 1, 2x) = r(l + u)e 



(f)~"W, 



U(v + 1/2, 2v + 1, 2x) = —= {2x)~ v K v {x) 

V7T 



(A4) 



where I v (x) and K v {x) are the modified Bessel functions. 
Then, taking into account that 



h/2{x) 



shr, K 1/2 (x) 



TVX ' V 2x 

and returning to the variable z, we obtain finally 

„ _ -Az -(z-z ) V \ 2 -a 



(A5) 



(A6) 



The latter equation describes the function u above the po- 
lar gap, z > zo, which results from continuous passage 
through the infinitesimal transverse current sheet accord- 
ing to Eq. ([6]). At z = zo we have u = e _Az °, and at z 3> zo 
Eq. (|A6|) takes the form 



(A7) 
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Pulsar force- free magneto sphere close to the magnetic axis 



log 4> 




log z 

Figure 1. The axial magnetic flux as function of the altitude above the neutron star based on Eq. J9J|. The approximations at z 
and z — y oo are shown by dots. 
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Figure 2. Field lines of the force-free magnetic field in the axial region of a pulsar (solid lines). Dotted lines correspond to the dipolar 
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